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Abstract 

In this paper, a generic self-similar flux rope model is proposed to 
probe the internal state of CMEs in order to understand the thermody- 
namic process and expansion of CMEs in interplanetary space. Using 
this model, three physical parameters and their variations with heliocen- 
tric distance can be inferred based on coronagraph observations of CMEs' 
propagation and expansion. One is the polytropic index V of the CME 
plasma, and the other two are the average Lorentz force and the thermal 
pressure force inside CMEs. By applying the model to the 2007 October 
8 CME observed by STEREO/SECCHI, we find that (1) the polytropic 
index of the CME plasma increased from initially 1.24 to more than 1.35 
quickly, and then slowly decreased to about 1.34; it suggests that there 
be continuously heat injected/converted into the CME plasma and the 
value of r tends to be |, a critical value inferred from the model for a 
force-free flux rope; (2) the Lorentz force directed inward while the ther- 
mal pressure force outward, and both of them decreased rapidly as the 
CME moved out; the direction of the two forces reveals that the thermal 
pressure force is the internal driver of the CME expansion whereas the 
Lorentz force prevented the CME from expanding. Some limitations of 
the model and approximations are discussed meanwhile. 



1 Introduction 

Coronal mass ejections (CMEs) are the most energetic eruptive phenomenon 
occurring in the Sun's atmosphere and the major driver of space weather. They 
carry a huge amount of mass, kinetic energy and magnetic flux into the inter- 
planetary space, and therefore may cause many significant consequences in the 
geospace. In this paper, we develop a generic flux rope model to infer three 
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physical parameters of CMEs as well as their variations with heliocentric dis- 
tance through the usage of the latest STEREO (Solar TErrestrial RElations 
Observatory) observations. The first parameter, is the polytropic index, T, that 
describes the thermodynamic process of the gas, and the other two are the 
Lorentz force and the thermal pressure force, that reveal the internal cause of 
the CME expansion. 

CMEs have been observed and studied for decades. There have been many 
observations, either through remote sensing observations or in-situ samplings, 
revealing the internal properties of CME plasmas. For example, the r emote sens 



ing da ta from SOHO/UVCS (ultraviolet coronagraph spectrometers. I Kohl et al 



2006() can diagnose the plasma temperature, density, velocity and heating at 
a few solar radii from the Sun. Such spectroscopic an alyses suggested that 
CMEs be a loop-like structure [e.g ., Ciaravella et aLl.l2003l with helical magnetic 



field [e.g., \Antonucci et all Il997t \Ciaravella et all l200fj, and probably have a 



higher temperature tha n that in the typical Corona in the near Sun region 



Ciaravella et al , 2003j . The thermal energy deposited into CME plasmas is 



roughly comparable to th e kinetic and gravitational potential energi es of CMEs 



in the inner corona [e.g., \Akmal et all 120011 : \Ciaravella et all |2001] . Some In- 



ternal properties of CMEs can also be revealed from in-situ observations, e.g, by 
Ulysses and ACE spacecraft. For example, the interplanetary CMEs (ICMEs) 
at 1 AU usually show a lower tem perature and stronger magnetic fields than 
that i n the ambient s o lar wi nd [e.g., Burlaaa et aUll981 ; I Richardson and CaneY 



19951: \Farruaia et all, 1 19931 . The ion char ge states in ICMEs are often higher 
[e.g., Lepri et all . 20011 Lynch et al. . 2003| . Based on the analysis of ion charge 
states, it was also found that the thermal energy i nput to the CME plas mas is 



Rakowski et all 12007 



at the same order of the CME kinetic energy [e.g. 

The above studies provided the information of the internal properties of 
CMEs, but only at a certain position and/or at a certain time. What is largely 
lacking is the global observations and thus the global understanding of the 
evolution of the internal state of CMEs during their continuous propagation 
throughout the interplanetary medium. What thermodynamic process does 
the CME plasma undergo? What happens with the various forces involved? 
Limited knowledge on these global issues were obtained through indirect ways, 
largely from the statistical combination of observations of many CMEs from 
multiple spacecraft over a long time period. The multiple-spacecraft measure- 
ments suggested that the polytropic index of CME plasmas be below 1.3 [e.g., 



Liu et all 12006 1 , which is appare ntly different from that of solar wind, which 



is about 1.46 Totten et al. I. ll995| . The radial widths of CMEs continuously in- 
crease at th e order of local Alfven speed a s they move away from the Sun within 
10 AU [e.g., I Wana and Richardsonl 12004 1 Wang et all 120051: \Jian et qZ.L |2008| . 
and the magnetic field decreases fa ster in ICMEs than in ambient solar wind 



but the temperature does not [e.g., I Wang and Richardson] . |2004 I Wana et al. 



20051 \Liu etal 



2006j . cd that all the above conclusions were 

established on statistical surveys, which can not review the detailed evolution 
behavior of any individual CMEs. 

It is now well accepted that CMEs, at least a significant percentage of them, 
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have a flux rope- like structure (Fig|T|) . Thus we try to study the internal state 
of CMEs by establishing a flux rope model. There are alrea dy various flux 



rope models concerning CME initiation and/or propagation [e.g., Burlaaa et al. 


1981 


; Goldstein, 1983; Cheri. 1989; Forbes and Isenbera. 1991; Kumar and Rust 


1996 


; Vandas et al.. 1997t Gibson and Loul 1998; Titov and Demoulin. 19991 



These models have their own specific purposes, and may not suit the issues 
attacked in this paper. We present our model in the next section. We then make 
a case study in Secf3] by applying the model to the CME that occurred on 2007 
October 8, whose expansion and propagation over a large distance throughout 
the interplanetary space were well observed. In Sec|H a brief summary is given. 
Finally, we thoroughly discuss the limitations and approximations of the model 
in Sec(5j 

2 The Model 

2.1 Derivation and Parameters 




Figure 1: Upper panel: Schematic plot of a flux-rope CME. Lower panel: The 
cross-section of the flux rope. R is the radius, L is the distance from the flux rope 
axis to the solar surface, and / is the axial length of the flux rope. The distance, 
h, of the CME leading edge and the minimum and the maximum position angles 
(indicated by the dot-dashed lines) are the three directly measurable quantities 
from coronagraph observations. 

As commonly assumed, CMEs are approximated as a cylindrical flux rope in 
a local scale, even though in the global scale they may be a loop-like structure 
with two ends rooted on the surface of the Sun. The flux rope can therefore be 
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treated axisymmetrically in the cylindrical coordinates (r, (j>, z) with the origin 
on the axis (ref to FigQ]) and ^ = J| = 0. As a result, only a 2-dimensional 
circular cross-section of the flux rope needs to be considered in our model. Let 
the radius of the circular cross-section be R, the expansion speed of the flux 
rope is given by 

,(.) - « a, 

Further, assuming that the flux rope is undergoing a self-similar expansion 
within the cross section, we can get a dimensionless variable 

T 

which is the normalized radial distance from the flux rope axis and independent 
of time. The x = 1 is the boundary of the flux rope. Therefore, the r-component 
of the velocity is 

v r (t,x) = ^ = xv e (t) (3) 

and the acceleration is 

orb*) = f^.* = ^(t)-%> (4) 
at xR 

where 

dv e (t) 

a e (t) = — (5) 

is the acceleration of the expansion. The first term on the right-hand side of 
EqfJ]is the acceleration of the radial motion of the plasma, and the second term 
is the acceleration contributed from the circular/poloidal motion. 

With the above preliminary preparation, we now investigate an arbitrary 
fluid clement in the flux rope by starting from the momentum conservation 
equation (in the frame frozen- in with the moving flux rope), 

^+p(v-V)v=-Vp+jxB (6) 

where p is the density, p is the plasma thermal pressure, B = (0, B,p,B z ) is 
the magnetic induction, and j = -^V x B is the current density. Here, the 
viscous stress tensor, S, gravity, F g , and the equivalent force due to the use of 
a non-inertial reference frame, F a opposite to the acceleration, are ignored (the 
validation of this treatment will be discussed at the end of this paper). This 
equation can be decomposed into the f and <f> components as follows 

dv r , dv r v, dp ,. _ s 
0: p_+p(^ + _i)=0 (8) 
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According to the self-similar assumption, v<j, has the following form 

v<f>(t,x) = f p (x)v p (t) (9) 

in which f p is a function of only x and v p is a function of only t. Combine it 
with Eq[5J and [5J it is inferred that 

v<t> = fc 1 x fe2 - 1 i?- fc2 (10) 

where k\ and ki are both constants. It is required that > 1 to guarantee 
that v<j, is physically meaningful, k% — 1 implies that the angular momentum 
of the flux rope is conserved, and ki > 1 means that the angular momentum 
decreases as the flux rope expands. Combine Eg El [51 [71 and [TUl we can rewrite 
the momentum conservation equation with f as 

(j x B) r = p(a e x - fc^-^-V 2 *'- 3 ) + R- 1 ^- (11) 

For a thermodynamic process, we relate the thermal pressure p with the density 
p by the polytropic equation of state 

V = k 3 p T (12) 

where k 3 is a positive constant and T is a variable treated as the polytropic 
index, and Eg [TT1 becomes 

(j x B) r = p(a e x - klR-^^x 2 ^-*) + hR- 1 ^- (13) 

ox 

Define a quantity f em to be the average Lorentz force over f from the axis to 
the boundary of the flux rope, f em — J Q (j x B) r c?r. From Eq[T51 we get 



fem ^< 



f pxdx-kjR- 2 ^- 1 f px 2ka - 3 dx + kaR- 1 f ^-dx (14) 
Jo Jo ® x 

fem > means that the average Lorentz force directs outward from the axis 
of the flux rope, causing expansion. On the other hand, f em < prevents the 
expansion of the flux rope. 

We assume that the mass of a CME is conserved when it propagates in the 
outer corona and interplanetary space, where the CME has fully developed. The 
mass conservation gives 

J prdrd(j)dz = 2nlR 2 J pxdx — M (constant) (15) 

where I is the axial length of the flux rope (FigQ]). Since the flux rope is assumed 
to be self-similar and it is generally accepted that the magnetic field lines are 
frozen-in with the plasma flows in corona/interplanetary space, the density in 
the flux rope has a fixed distribution f p (x), and therefore 

p{t,x) = f p {x) Po {t) (16) 
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Define positive constants 



k 4 = / fpxdx (17) 
Jo 

k 5 = [ f p x 2k *- 3 dx (18) 
Jo 

and a variable 

9(r) = / p r (o)-/ p r (i) (19) 

Then it can be inferred from Eq[T5lthat 

Po = ^-k^MR-H- 1 (20) 

Z7T 

and EqfTHcan be written as 

fern = ^K^ 2 /- 1 - k\k^ fc 5 iT 2fe2 ~H~ 1 ) - f th (21) 
Z7T 



fth = i f R -^dr = k 3 qp T Q R- 1 (22) 



where 



is the average thermal pressure force. Like f em , fth points outward if it is larger 
than zero. 

On the other hand, in an axisymmetric cylindrical flux rope, 

B = B^cf) + B z z = VxA (23) 

B* = ~ d -t ( 24 ) 
B z = i&M,) (25) 

As the magnetic flux is conserved in both <p and z directions, we get 

n = -lf R§ tdr = l(A z (0)-MR)) (26) 
= 2tt f Q R ^{rA^dr = 2nRA^R) (27) 

In order to satisfy the self-similar expansion assumption, A^ and A z have to 
keep their own distributions, respectively. Thus, according to the above two 
equations, 

A,(t,x) = (28) 
A z (t,x) = ^ (29) 
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It can be proved that the conservation of hclicity is satisfied automatically 

H m = J B • Ardrd<j)dz 

f 1 f d df 

= constant (30) 
03 [2H1 and E51 we can calculate the Lorentz force in the flux 



Combining Eqfj 
rope 

jxB 



— (V x B) x B 

Mo 



dx dx dx 



(31) 



and therefore 

fern — 



-Mo 1 



^ R 1 h x -dx"d-x {x -dx- )dx 



-fi 1 k 6 R~ 



(32) 



where and k-? are both constants. It could be proved that the sign of ke is 
determined by B 2 (R) - B 2 (0), and k 7 > 0. 

The two forms of f em , Eq|2T] and [32j result in 

= M- x (k 6 R-H + k-iR-H- 1 ) + 2irM- 1 k 3 (2irk 4 M- 1 R 2 l)~ r [fi (0) - / p r (l)]7fS3) 

in which f t h is substituted by Eq |22l As at present it is impossible to practically 
detect the axial length of a flux rope, here we will relate it with a measurable 
variable, L, the distance between the flux rope axis and the solar surface (Fig{T|) , 
at which altitude the flux rope originates, by the assumption 



/ = kxL 



(34) 



where k$ is a positive constant. The topology of flux rope as shown in Figure [T] 
implies that this assumption is reasonable. Finally, Eg 1331 can be simplified to 



where 



c R' 



= -c 2 R- 3 L - c.R-'L- 1 + c 4 (c r 5 - cDR^L 1 - 7 
c = fc^fc^ 1 fc 5 > 



(35) 
(36) 
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Table 1: List of the constants k%s and cq-6 



Constant Interpretation 


Constant Interpretation 


k\ Scale the initial magnitude 
of the poloidal motion 

&2 Decrease rate of the angular 
momentum as the flux rope 
expands 

fc 3 Coefficient in the polytropic 
equation of state 


hi and k$ Integral constants related 
to the density distribution 

ks and k-j Scale the initial magnitude 
of the Lorentz force con- 
tributed by the axial and 
poloidal fields 
k$ Assumed constant to relate 
the length of flux rope I to 
distance L 


Co Scale the initial magnitude 
of the acceleration due to 
the poloidal motion 

ci Similar to ki 


C2 and C3 Similar to fcg and k 7 

C4(c5 — Cg) Scale the initial magnitude 
of the contribution by ther- 
mal pressure force 



ci = 2k 2 - 2 > (37) 

c 2 = 2ir^ 1 M- 1 k 6 k 8 (38) 

c 3 = 2irnQ 1 M- 1 k 7 kg 1 > (39) 

c 4 = 2TrM- 1 k 3 k s > (40) 

c 5 = (2^)- 1 Affc 4 - 1 fc 8 - 1 / P (0) > (41) 

c 6 = (2n)- 1 Mk^ 1 k 8 1 f p (l) > (42) 

The left-hand side of Eql35] describes the motion of the fluids in the flux rope. 
Its first item is the acceleration due to the radial motion (i.e., expansion) and 
the second one gives the acceleration due to the poloidal motion. The right-hand 
side reflects the contributions from the Lorentz force (the first two items) and 
thermal pressure force (the last one). The constants an d c o-6 appeared 
above are summarized in Table [1] 

The Lorentz force and thermal pressure force can be rewritten in terms of 
the constants cq-q, fc§ and the total mass M as follows 

M K , o 

fern = ~7T^^ R + L ) (43) 



27rfc 8 
M 
2Trk R 



f* = -£r*(£-$)Br ,a '- 1 L-r (44) 



and their ratio is 



f ^ »2r-4rr 1 _ n2r-2rr-2 

Jem C 2 H L + C 3 K L 

-f— = 7T — n ( 45 ) 

fth c 4 (c 6 -c^) 

In summary, starting from MHD equations with the three major assumptions 
that (1) the flux-rope CME has an axisymmetric cylinder configuration, (2) its 
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cross-section is self-similarly evolving, and (3) its axial length is proportional to 
the distance from the solar surface, we find that the polytropic index, T, can 
be related to the measurable parameters: the distance, L, the radius, R, and 
another derived quantity, the expansion acceleration (a e ), as shown in Eq |35l If 
we have enough measurement points, the unknown constants co-6 an d variable 
T could be obtained through some fitting techniques (e.g., that described in the 
first paragraph of Sec l3.2p , and then the relative strength of the Lorentz force 
and thermal pressure force can also be easily calculated by Eq03] and 2U 



2.2 Asymptotic Value of Polytropic Index T 

Here, we consider the case of a nearly force-free expanding flux rope. It is 
generally true that most CMEs are almost force-free at least near 1 AU though 
they may be far away from a froce-free state at initial stage. It can be proved 
that R oc L (ref. to Appendix), i.e., R = aL where a is a positive constant. 
Then Eq |4"5"1 becomes 

f (n ™ 2r - 4 j. n ™ 2T - 2 \ r3r-4 
Jem _ (,c 2 a + c 3 a )L 

fth c 4 {c L 6 -c L 5 ) 

It is found that r = | is a critical point, above/below which the absolute 
value of Lorentz force decreases slower/faster than that of thermal pressure 
force as inc reasi ng distance L. This valu e of T is the same as that obtained by 
Zowl [19821 ] and \Kumar and RusA |l996l | for a self-similar expanding flux rope. 



This inference is reasonable because smaller T implies the plasma absorbs more 
heat for the same expansion and therefore the thermal pressure should decrease 
slower. 

Under force-free condition, Eq[35]can also reduce to 

a e = c„{aLy ci - 3 - c 2 a- 3 L- 2 - c 3 a^L- 2 + c 4 (c£ - c^a 1 -^ L 2 ~ 3r (47) 

and at infinite distance, L — > +oo, we have 

a eoo - - cl^a 1 -^-^ 3r ~ (48) 

The above equation indicates that r = | is another critical point. The poly- 
tropic index T should be larger than | to make sure that the flux rope will 
finally approach a steady expansion and propagation state (including the case 
that the flux rope stop somewhere without expansion). Otherwise, the flux rope 
will always accelerated expanding. 

Based on the current observations, the expansion behavior of CMEs at large 
heliocentric distance is not as clear as that in the inner heliosphcrc. The in- 
vestigations on the radial widths o f CMEs suggest that CMEs at least keep 



expan ding within about 15 AU [e.s. AWano and Richardsom . 12004 1 Wang et al. 



2005l | . but the expansion speeds seem to be slower and slower. Although the 



number of CMEs observed near and beyond 15 AU is small and the uncertainty 
of statistics is large, it is likely that a CME may not be able to keep an accel- 
erated expansion always. Thus, in practice, the polytropic index of the CME 
plasma should be larger than 2. 
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3 The 2007 October 8 CME 

3.1 Observations and Measurements 

The suite of SECCHI instruments on board STEREO spacecraft provide an 
unprecedented continuous view of CMEs from the surface of the Sun through 
the inner heliosphere. The instruments, EUVI, COR1, COR2, HI1 and HI2, 
make the images of the solar corona in the ranges of 0-1.5, 1.4-4 . 0, 2.5 -15.0, 
-15-90, and ~90-300 solar radii (R s ), respectively [Howard et all l2008j |. The 



SECCHI observations present us the great opportunity to study the evolution 
of CMEs over an extended distance. The CME launched on 2007 October 8 
is a well observed event, which is used to study the CME evolution and the 
applicability of our model. 

This CME was initiated close to the western limb as seen from STEREO 
B. Hereafter all the observations used are from instruments on board the B 
spacecraft. Figure [2] shows five images of the CME at different distances from 
the Sun. The CME was accompanied by a prominence eruption starting at 
about 07:00 UT on October 8, as seen by EUVI. The CME source region is 
clearly shown in the EUVI 304A image on the top-left panel of the hgure. The 
erupting prominence was also seen in the COR1 running-difference image (the 
top-right panel). The CME was first observed in COR1 at about 08:46 UT on 
October 8, and continuously ran through COR2 and HI1 fields of view (FOV). 
It even showed in the HI2 FOV after about 12:00 UT on October 10. Since the 
CME was launched from the western limb and showed a circular-like structure, 
we believe that the CME was viewed by the instruments through an axial-view 
angle. Therefore, the projection of the CME on the plane of the sky can be 
treated as the cross-section of the CME. 

To obtain the two quantities, R and L, for necessary model inputs, here we 
simply measure three parameters, the heliocentric distance of the CME leading 
edge, h, and the maximum and minimum position angles, 9 ma x and 9 m im of the 
CME as shown in Figure [1] Under the assumption of a circular cross-section, 
R and L could be derived by 

R = h-L-Rs (49) 

L = -l , • oL ~ R s (5°) 
1 + sin '" ax 2 """ 

It should be noted that the measurements in HI2 images arc not included in the 
following analysis, because the elongation effect is not negligible. 

Figure [3] shows the measurements and the derived parameters. The CME is 
a slow and gradually accelerated event. It took about 46 hours for its leading 
edge to reach 70 Rs- Nevertheless, because of its slowness, we are able to make 
about one hundred measurement points for this CME. The red crosses plotted in 
the leftupper panel suggests that the CME angular width increased at the early 
phase (mainly in the COR1 FOV), and then reached to a near-constant value in 
the COR2 and HI1 FOVs. The right panel presents the evolution of the derived 
R and L. It is shown that the radius of the flux-rope CME is about 20 Rs when 
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Figure 2: Images of the 2007 October 8 CME taken by (a) EUVI at 304A, 
(b) COR1, (c) COR2, (d) HI1 and (e) HI2 on board STEREO B. The panel at 
lower-right corner shows a direct image of the CME in the HI1 FOV. This image 
has been corrected to the plane perpendicular to the line between the Sun and 
STEREO B, because the CME is assumed to be a limb event in CORs and the 
direction that HI1 camera faces to is different with CORs'. 
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16:00 00:00 08:00 16:00 00:00 06:00 



Figure 3: Left-upper panel: The measurements of the heliocentric distance, h, of 
the flux-rope CME leading edge and its angular width, A9 = max — 8 m i n . Right 
panel: The derived distance, L, of the flux rope axis from the solar surface and 
the flux rope radius, R. Left-lower panel: The propagation, v c , and expansion, 
v e , speeds derived from L and R, respectively. 



it propagated nearly 50 Rs away from the Sun, which put the leading edge at 
about 70 Rs- The left-lower panel exhibits the speeds derived from the R and 
L, namely expansion, v e and propagation, v c , speeds, respectively. At the early 
phase, the expansion speed was very close to the propagation speed. In the 
later phase, the propagation speed increased more quickly than the expansion 
speed. The increased difference between v c and v e is probably because of (1) 
the enhanced drag force of the ambient solar wind, which is fully formed in the 
outer corona and (2) the weakened pressure in the CME. The issue of CME 
acceleration, which is as important as CME expansion, is not addressed in our 
model presented in this paper. 

In the measurements, the CME radius obtained is the one along the latitu- 
dinal direction on the meridional plane. This radius would be the same as the 
radius along the radial direction if the cross-section is a perfect circle. However, 
the true cross-section deviates from the perfect ci rcle, and the deviat i on be comes 



larger as the CME is further from the Sun [e.g.. li?^e^ and Crooken . 120041 ] . The 
distortional stretching of the cross-section is caused by the divergent radial ex- 
pansion of the background solar wind, which causes kinematic expansion of 
CMEs along both the meridional and azimuthal directions, but not at all along 
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the radial direction. The CME expansion along the radial direction is mostly 
driven by the dynamic effect, such as pressure gradient forces, while the expan- 
sion along the other two directions that lie on the spherical surface is caused by 
the combination of the dynamic and kinematic effec ts. As a result, the overal l 



cross-section is a convex-outward "pancake" shape [Rilev and Crooken . 12004 1 . 
Figure [2f shows such a distortion of the 2008 November 8 CME as observed in 
HI1 FOV; the aspect ratio, defined by the ratio of the radius along the merid- 
ional direction and that along the radial direction, is about 1.4 when the CME 
leading edge is at ~ 70Rs- 

Due to this stretching effect, our measurements assuming a circular cross- 
section lead to the inaccuracy of the measured parameters and the inferred 
parameters as well. In order to study the internal state of a CME, the radius 
of the CME, R, should be the one along the radial direction, and it is appar- 
ently overestimated when the radius along the meridional direction is adopted. 
The derived expansion speed of CME is thus larger than the true value. Such 
simplified measurements would infer unrealistic parameters of CME at 1 AU. 
For instance, the observed radius of 20 Rs of the CME at a distance of 50 Rs 
from the Sun would imply a CME cross-section of 0.8 AU at 1 AU, which is 
too larger to be true. The observed speeds of v e and v c would imply a speed of 
about 150 km/s at the trailing edge of the CME, which is much smaller than 
the observed solar wind speed, i.e., about 300 km/s. Therefore, one should be 
cautious when our method is applied to CMEs at a large distance from the Sun 
(e.g, > 70i?s). The heliospheric region we investigated in this paper is within 
about 70 Rs, and the stretching effect is relatively small. Nevertheless, we will 
carefully estimate the errors on CME parameters in the second paragraph of 
SecJS] We point out here that there is an observational difficulty in measuring 
the radius of CMEs along the radial direction in a consistent way, mainly be- 
cause of the low brightness contrast of the CME trailing edge in coronagraph 
images. This difficulty might be overcome if the CME of interest is particularly 
bright. 

Before modeling the CME, we fit the measurement points with a certain 
function to retrieve the smooth evolution process of the CME, which is required 
for the model. We use the modified function of log normal distribution to fit 
the speeds. We did not fit the expansion acceleration directly, because any 
small error in measurements of R will be dramatically amplified in its second 
derivative a e . The fitting function of velocity has the form 



v(t) = — 



1 + erf 



ln(i-fo) -M 



SV2 

where erf (z) is the erf function or error function, defined by 



(51) 



erf(z) ee I e't'dt (52) 



in 



This function has a value range from to . It is chosen because the measure- 
ments show a trend that, at least within the FOVs of SECCHI, both the speeds 
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will not increase forever, but instead asymptotically approach a constant speed, 
Woo ■ The acceleration can be derived by 

, , Woo [ln(t-tp)-M] 2 

a e (t) = — =2S e ^ 53 

The solid lines in the left-lower panel of Figure [3] show the fitting results. The 
fitted parameter is 118 km/s for expansion and 246 km/s for propagation. 
As a comparison with the measurements, the integrals of the fitting curves of 
the speeds are also plotted in the right panel. It has been mentioned before that 
these estimated speeds suffer the solar wind stretching effect. Particularly, the 
estimated expansion speed is larger than that it should be. The error will be 
discussed in SecEl 

3.2 Results 

To fit the above curves with the model, Eq |35[ we use an iterative method. Gen- 
erally speaking, first we solve this equation in every 8 neighboring measurement 
points to obtain a set of parameters cq-q and T. The segment of the 8 points 
is a running box through the entire evolution process of the CME. Secondly, 
input the obtained variable T into the model as guess values to fit the global 
constants co-6- Thirdly, use the fitted co-6 to update the variable T by solving 
Eq[35] again. Then iterate the above 2nd and 3rd steps to make constants co-6 
and r converging to a steady solution. For the sake of simplicity, we ignore the 
poloidal motion of the fluid by setting Co zero. It is also because there seems no 
strong observational evidence showing a ring flow inside a CME. 

The model results are shown in Figure |4] The uncertainty of the model 
results is estimated from the relative error of a e , which is given by 



E = 



(54) 



where a e i is the value calculated by the input data, and a em is the model value. 
The error curve is plotted in the left-upper panel of Figure [4j It is found 
that the error is smaller than 1%, except during 12:00 - 18:00 UT. A possible 
explanation of the large uncertainty during that time has been given in the last 
second paragraph of Sec [5] 

3.2.1 Polytropic Index 

From the right- upper panel of Figured! it is found that T was less than 1.4 
throughout the interplanetary space. In the inner corona, say L < 2i?s, it 
was about 1.24. After entering the outer corona, it quickly increased to above 
1.35 at L ss 5i?s, and then slowly approached down to about 1.336, which is 
very close to the first critical poi nt |. Thi s value of T is consistent with the 
observational value obtained from Liu et al. 2006j statistics for protons. As the 



CME kept expanding during its propagation in the FOVs, the polytropic index 
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less than | means that there must be some mechanisms to inject heat from 
somewhere into the CME plasma. Although the CME plasma continuously got 
thermal energy, the proton temperature may be still much lower than that in 
the ambient solar w ind, as revealed by many in-situ observations of MCs [e.g., 
Burlaaa et all Il981 |. 




16:00 00:00 08:00 16:00 00:00 08:00 
Start Time ( B-Oct-2007 09:00:00) 



00 00:00 08:00 16:00 00:00 08:00 
Start Time ( 8-Oct-2007 09:00:00) 



Figure 4: Left-upper panel: The profile of a e (black), the modeled result (dashed 
green), and the relative error (see text for details). Right-upper panel: The 
variation of the polytropic index, T, of the CME plasma. Left-lower panel: The 
variations of the average Lorentz force, f em , and the average thermal pressure 
force, fth- Their signs have been marked on the upper right corner. Right-lower 
panel: The sum and ratio of the two forces. 



We believe that the hot plasmas in the lower solar atmosphere is probably a 
major heat source of CMEs in the interplanetary space. As shown in Figure [TJ 
a CME is believed to be a looped structure with two ends rooted on the solar 
surface i n a global scale . Bidirectional electron str eams are one of the evidence 
of it [e.g. Farruaia et al. . 1993I : Larson et al . 1997 . Thus it is possible that heat 
is conducted from the bottom to CMEs. The ambient solar wind with higher 
temperature might be an additional source because the temperature difference 
between the two mediums is significant. However, the cross-field diffusion of par- 
ticles are much more difficult than the motion parallel to magnetic field lines, 
especially in a nearly force- free flux rope; the coefficient ratio k±/k\i of perpen- 
dicular to parallel diffusion roughly locates in the range of 0.005 — 0.05 [e.g. , 



Jokiv ii et all\l99fx\Giacalone and Jokivin . 1999; Zank et all\200<x\Bieber et al. 



2004j . Thus the contribution of the ambient high-temperature solar wind should 



be very limited. 

It is well known that the magnetic energy decreases as CMEs propagate 
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away from the Sun. According to our model, the total magnetic energy is given 

by 



m 2 M o / 




B 2 rdrd(j)dz = nfi 1 {k$l~ 1 + k w R~ 2 l) 



(55) 



where 



(56) 
(57) 



are both positive integral constants. The magnetic energy generally dissipates 
at the rate of ~ which is a significant dissipation as CMEs move outward. 
However, such magnetic energy dissipation does not necessarily mean to be a 
major source of the heat. According to MHD theory, magnetic energy partially 
goes into kinetic energy and partially converts to thermal energy. The former is 
due to the work done by Lorentz force (j x B • u) , and the latter is through the 

-2 

Joule heating (^-) process, where j is the current density and a is the electrical 
conductivity. Since a usually has a large value in interplanetary medium, with- 
out anomalous resistivity, the magnetic energy does not have an efficient way 
to be converted to thermal energy. However, there are possibly many non-ideal 
processes, such as turbulence, but not accounted by MHD theory. Thus we do 
not know whether the dissipated magnetic energy is a major source of heating 
or not. 

3.2.2 Internal Forces 

The averaged Lorentz force, f em , and thermal pressure force, fth, have been 
presented in the left-lower panel of Figure [4j Their absolute values are very 
close to each other, and both of them decreased continuously throughout the 
interplanetary space. The signs of the two forces are opposite. f em is negative 
indicating a centripetal force, whereas fth is positive, centrifugal. This result 
suggests that the thermal pressure force contributed to the CME expansion, but 
the Lorentz force prevented the CME from expanding. 

The difference between the two forces can be seen more clear from the right- 
lower panel of Figured] The black line exhibits the net force, fem+fth, inside the 
CME. ft directed outward and reached the maximum at about 10:30 UT. The 
profile is consistent with the expansion acceleration presented in the left-upper 
panel (the black line). Thus the net force just shows us the internal cause of 
the CME expansion. The red line is the ratio of their absolute values. Its value 
changed in a very small range from about 1.0 to 0.98. It suggests that such a 
small difference between the two forces is able to drive the CME expanding with 
the acceleration at the order of 1 m/s 2 . Moreover, the ratio decrease means that 
the Lorentz force decreased slightly faster than the thermal pressure force. One 
may notice that, since V was larger than the first critical point | at L > 5i?s, 
according to the analysis in Sec l2.2[ the Lorentz force should drop slower than 
thermal pressure force. Actually it may not be an inconsistency, because the 
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inference derived in Sec l2.2l is established on the force- free assumption, the CME 
we studied may not be perfectly force-free, and therefore the first critical point 
of r probably shifts a little bit. 

Usually, CMEs are a flux rope with two ends rooted on the Sun. The axial 
curvature of the flux rope may cause the magnetic strength at the Sun-side of the 
flux rope larger than that at the opposite side, which leads the Lorentz force hav- 
ing an additi onal component to drive the flux rope moving outward away from 



the Sun \e.s..\Garren and Jamed . \1994k\Lin et aZ.l .ll998. 2002; Kl iem and Torok 



20061: [Fan and Thus, as the flux rope we applied here is assumed 



to be a straight cylinder, the Lorentz force f em we derived does not include 
the component caused by the axial curvature of the flux rope. This component 
is important in studying the propagation properties of a CME. However, our 
model is to study the CME internal state (specifically the thermodynamic pro- 
cess and expansion behavior), and its propagation behavior is obtained directly 
from coronagraph observations, thus the neglect of this component should be ac- 
ceptable although it does bring on some error, which has been briefly mentioned 
in the second paragraph of Sec 03 



4 Summary 

In this paper, we developed an analytical flux rope model for the purpose of 
probing the internal state of CMEs and understanding its expansion behavior. 
The model suggests that, if the flux rope is force free, there are two critical values 
for the polytropic index T. One is |, above/below which the absolute value 
of the Lorentz force decreases slower/faster than that of the thermal pressure 
force as the flux-rope CME propagates away from the Sun. The other is |, 
above which the flux-rope CME will essentially approach a steady expansion 
and propagation state. 

By applying this model to the 2007 October 8 CME event, we find that (1) 
the polytropic index of the CME plasma increased from initially 1.24 to more 
than 1.35 quickly, and then slowly decreased to about 1.336; it suggests that 
there be continuously heat injected/converted into the CME plasma and the 
value of T tends to be the first critical value |; (2) the Lorentz force directed 
inward while the thermal pressure force outward, both of them decreased rapidly 
as the CME moved out, and the small difference between them is consistent with 
the expansion acceleration of the CME; the direction of the two forces reveal 
that the thermal pressure force is the internal driver of the CME expansion, 
whereas the Lorentz force prevented the CME from expanding. 



5 Discussion 

In our model, the interaction between CMEs and the solar wind has been implic- 
itly included to certain extent, though we do not explicitly address these effects. 
The consequences of the interaction, in terms of the effects on the CME dynamic 
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evolution can be roughly classified into the following three types: (1) the solar 
wind dragging effect, which is is due to the momentum exchange between the 
CME plasma and the ambient solar wind and mainly affects the CME's propa- 
gation speed or the bulk motion speed, (2) the solar wind constraint effect on 
expansion, which is caused by the presence of the external magnetic and thermal 
pressures and mainly prevents a free expansion of the CME (i.e., in all direc- 
tions), and (3) the solar wind stretching effect on expansion, which is caused 
by the divergent radial expansion of solar wind flow, and causes flattening or 
"pancaking" of CMEs. The first two effects are indirectly included in the model 
through the measurements of L and R. Different dragging and/or constraint 
force(s) may result in different variation of L and/or R with time (or heliocentric 
distance) . Particularly, we do not need to explicitly put the solar wind dragging 
term in the model, because we are addressing the internal state of CMEs, not 
the bulk acceleration. The stretching effect, which is of a kinematic effect, is 
not included in our model. As discussed in the sixth paragraph of Sec J3.ll this 
is largely due to the limitation of the measurements. The possible errors caused 
by such effect are explicitly addressed in the next paragraph. 

The main uncertainty of this model, we believe, comes from the assump- 
tion of an axisymmetric cylinder, in which the curvature of the axis of the 
flux rope and the distortion of the circular cross-section are not taken into ac- 
count. As to the first one, the neglect of the axial curvature generally results 
in the Lorentz force underestimated. As to the second one, as discussed ear- 
lier, the distortion of the CME cross-section is due t o the kinematic stretching 



effect of a spherically divergent solar wind flow [e.g. , Crooker and Intriliaator 
19961: Russell and Mulligan . 2002 : Riley et al . 2003 : \Rileu and CrookeA . Is~ 



2004 



Liu et ali . 1.2006] . In the case of the particular CME studied in this paper, the 
aspect ratio is about 1.4 when the CME leading edge is at ~ 70Rs (or the flux 
rope axis is at ~ 56i?s). The overall shape of the CME looks like an ellipse. 
To estimate the errors caused by the circular assumption, we approximate the 
ellipse to be a circle of the same area. With this treatment, we estimate that 
R is overestimated by 19%, and L is underestimated by 11%. Therefore, the 
expansion speed is overestimated by 19%, and the propagation speed is under- 
estimated by 11%. Further, we find that the density is underestimated by 21% 
(ref to Eq |20|) . f t h is underestimated by 39% (ref to Eg l22l and assume r = |), 
f em is underestimated by 25-58% (ref to Eq |3"2")l . and the error of the polytropic 
index is probably neglected (ref. to EqfJ3|). These errors are evaluated for the 
CME at ~ 70Rs- At a smaller distance, we expect that the errors be smaller, 
since the distortion is less severe. 

The self-similar assumption made in our model may be another error source, 
in which we assume that the distributions of the quantities along f in the flux 
rope remain unchanged during the CME propagates away from th e Sun . Self- 
similar evolution i s a frequently used assumption in modeling [e.g. , 



Kumar and \Gibson and Lowl.ll998t I Krall and St. Cvn , 



Low, 1982 



2006]. The 



recent research bv lDemoulin and Dasso\ [2009j suggested that, when I, the length 
of flux rope, is proportional to p t ^ 4 , the total pressure in the ambient solar 
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wind, a force-free flux rope will evolve self-similarly. The total pressure of solar 
wind consists of thermal pressure pth = nkT and magnetic pressure p m = . 
Near the Sun, we can assume that the magnetic pressure is dominant, thus it 



oc L. Since the length of 

-1/4 



is approximated that p t ~ p rn oc L 4 , i.e., p t 

a flux rope is usually proportional to the distance L, we have I oc p t ±, ' ± . It 
means that self-similar assumption should be a good approximation when the 
CME is nearly force-free and not too far away from the Sun. Other previous 
studies also showed that t he self-similar evo l ution of CMEs is probably true 



withi n tens solar radii [e.z. AChen et all Il997t \Krall et all 120011 : \Maricic et al. 



2004]. On the other hand, however, the self-similar assumption must be broken 
gradually. An obvious evidence is from the solar wind stretching effect as have 
been addressed before. Another evidence is that a CME may relax from a 
complex stru cture to a nearly fo rce-free flux rope structure, for example the 
simulation bv lLvnch eiai\ [2004j . 

For the CME plasma, neglecting the viscous stress tensor in Eq[6] might be 
appropriate. The viscous stress tensor of protons can be approximately given 
by 



BiBj 



B • B Vv V • v 



B 2 



(58) 



and 770 is th e coefficient of viscosity that cou ld be estimated by 770 « 10 T ; 



kg-m 1 -s 1 Braainskii . 1965 : Hollwea , 19851 ]. Here Sij is the unit tensor, v is 



flow velocity, and T p is the proton temperature. Since the proton temperature 
in CMEs is low, r/o and therefore the viscous stress tensor is very small. Thus, 
we guess that the viscosity in the momentum equation might be ignored. 

Both forces ignored in Eq[51 the gravity F s and the equivalent fictitious force 
F Q due to the use of a non-inertial reference frame, are in the radial direction 
in the solar frame. Their effects can be evaluated by comparing them with 
the acceleration of the expansion of the fluids in flux rope CMEs. The solar 
gravity acceleration is about 270 m/s 2 at the surface, and decreases at the rate 



of 



■, which makes it as low as ~2.7 m/s at 10 Rs- And F a should be also 



very small for most CMEs beyond 10 Rs- Thus both forces would significantly 
distort the model results only on CMEs with slow expansion acceleration in the 
lower corona, but not on those with large expansion acceleration or in the outer 
corona. This may be the reason why a large error of a e appears during 12:00 - 
18:00 UT in modeling this CME (left-upper panel of Figg}. 

The flux-rope model presented in this paper might be the first of its kind to 
provide a way to infer the inter state of CMEs directly based on coronagraph 
observation s. It is different from other C ME dynamic models, such as those by 
Chen 19891 ] and Gibson and Low 19981 ] , which were designed to study the in- 
teraction of CMEs with the ambient solar wind and other dynamic processes by 
adjusting the initial condi tions of CMEs and the global parameters of the am- 
bient solar wind. Besides, Kumar and Rust 19961 ] proposed a current-core flux 
rope model with self-similar evolution (ref. to KR model thereafter). Although 
a self-similar flux rope is also employed in their model, our model is largely 
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different from theirs. First , the flux rope in KR model is assumed force-free 
and the Lundquist solution Lundauist , 1950j ] is applied to describe the internal 
magnetic structure, but our model does not specify the magnetic field distri- 
bution and it may be non- force- free. Secondly, the self-similar assumption in 
KR model limits the radius of the flux rope to be proportional to the distance, 
whereas our self-similar condition is held only in the cross-section of the flux 
rope; the R and L in our model are two independent measurements (see Fig|3]). 
Thirdly, KR model did not consider the solar wind effects on the flux rope, while 
two of three solar wind effects are implicitly included in our model. Thus one 
can treat our model a more generic one. Undoubtedly, KR model is an excellent 
model for force-free flux ropes, and got many interesting results. For example, 
it is suggested that the polytropic index is | for a CME far from the Sun. It is 
an inference from their self-similar assumption, and it seems to be true for the 
2007 October CME we studied here. In our model, the V value of | implies a 
special case (Sec l2.2p in which the two internal forces f em and fth vary at the 
same rate. Further work will be performed to test whether it holds for all CME 
events. 
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Appendix 

In cylindrical coordi nate system, the magnetic field of a force-free flux rope has 
the Lundauist 1950j solution 



B r = 

B^ = H B Q J x (2.41a;) (59) 
B z = B J (2Alx) 

where x = is the normalized radial distance as defined in Sec[H Jo and J\ 
are the zero and first order Bessel functions, H = ±1 indicates the sign of the 
handedness and B$ is the magnetic field magnitude at the axis of the flux rope. 
According to the properties of Bessel function, we have the magnetic vector 
potential 

A f = ^ M2A1X) (60) 
A z = ^g^J (2.41x) (61) 
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The conservation of Q z 

<P z =2tt / — (rA^dr = Ji(2.41) = constant (62) 

Jo o r 2.41 

requires that 

B a = 2A\ ai R- 2 (63) 
where ai is a constant. The magnetic vector potential can be rewritten as 

A = ^Ji (2.41a;) (64) 
rt 

A z = ^J (2.41z) (65) 
K 

Meanwhile, the magnetic helicity is 

H m = J B • Adrrd<j)dz = A.Wl-Ka\o2HR~H (66) 

where 02 = f x( Jq + J^)dx is a constant. The conservation of if m results in 

i? oc I (67) 

Combined it with the assumption Eq |341 it is inferred that 

R oc L (68) 

which means that the force-free flux rope expands radially. 
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